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ABSTRACT 

We measure the distribution of carbon in the intergalactic medium as a function of redshift z and overdensity 5. 
Using a hydrodynamical simulation to link the H I absorption to the density and temperature of the absorbing 
gas, and a model for the UV background radiation, we convert ratios of C IV to H I pixel optical depths into 
carbon abundances. For the median metallicity this technique was described and tested in Paper I of this series. 
Here we generalize it to reconstruct the full probability distribution of the carbon abundance and apply it to 19 
high-quality quasar absorption spectra. We find that the carbon abundance is spatially highly inhomogeneous 
and is well-described by a lognormal distribution for fixed 6 and z. Using data in the range log 6 = -0.5 — 1.8 
and z = 1.8-4.1, and a renormalized version of the Haardt & Madau (2001) model for the UV background 
radiation from galaxies and quasars, we measure a median metallicity of [C/H] = -3.47^+0.08^(z-3) + 
0.65^];!2(log£-0.5) and a lognormal scatter of <r([C/H]) = 0.76^g;gf + 0.02^;?|(z-3)-0.23^(log^-0.5). 
Thus, we find significant trends with overdensity, but no evidence for evolution. These measurements imply that 
gas in this density range accounts for acosmic carbon abundance of [C/H] = -2.80 ±0.1 3 (ftc ~ 2 x 10~ 7 ), with 
no evidence for evolution. The dominant source of systematic error is the spectral shape of the UV background, 
with harder spectra yielding higher carbon abundances. While the systematic errors due to uncertainties in the 
spectral hardness may exceed the quoted statistical errors for 6 < 10, we stress that UV backgrounds that 
differ significantly from our fiducial model give unphysical results. The measured lognormal scatter is strictly 
independent of the spectral shape, provided the background radiation is uniform. We also present measurements 
of the C III/C IV ratio (which rule out temperatures high enough for collisional ionization to be important for 
the observed C IV) and of the evolution of the effective Lya optical depth. 

Subject headings: cosmology: miscellaneous — galaxies: formation — intergalactic medium — quasars: ab- 
sorption lines 



1. INTRODUCTION 

The enrichment of the intergalactic medium (IGM) with 
heavy elements provides us with a fossil record of past star 
formation, as well as a unique laboratory to study the effects 
of galactic winds and early generations of stars. With the 
advent of the HIRES spectrograph on the Keck telescope it 
rapidly became clear that a substantial fraction of the high col- 
umn density (Nm > 10 14 5 cm" 2 ) Lya absorption lines seen at 
redshift z ~ 3 in the spectra of distant quasars have detectable 
associated absorption by C IV (Cowie et al. 1995; Ellison et 
al. 2000) and Si IV (Songaila & Cowie 1996). More recently, 
observations with the UVES instrument on the Very Large 
Telescope (VLT), as well as with the Keck Telescope, have 
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revealed that the same is true for O VI at z ~ 2 (Carswell, 
Schaye, & Kim 2002; Simcoe et al. 2002; Bergeron et al. 
2002, see also Telfer et al. 2002). Both simple photoioniza- 
tion models (e.g., Cowie et al. 1995; Songaila & Cowie 1996; 
Carswell et al. 2002; Bergeron et al. 2002) and numerical sim- 
ulations assuming a uniform metallicity (Haehnelt et al. 1996; 
Rauch et al. 1997a; Hellsten et al. 1997; Dave et al. 1998) in- 
dicate that the gas has a typical metallicity of 10~ 3 to 10~ 2 
solar. 

The gas in which metals can be detected directly 
is thought to be significantly overdense [S = pj (p) ~ 
10 (Vm/10 15 cm" 2 ) 2 / 3 [(l+z)/4]- 3 ; Schaye 2001] and conse- 
quently fills only a small fraction of the volume. The typical 
metallicity of the low-density IGM (away from local sources 
of metals) remains largely unknown, although statistical anal- 
yses based on pixel optical depths do indicate that there is 
C IV (Cowie & Songaila 1998; Ellison et al. 2000) and O VI 
(Schaye et al. 2000a) associated with the low-column density 
Lya forest. 

Despite the recent progress, many questions remain regard- 
ing the distribution of metals in the IGM: Are metals generally 
present at S <C 10? Does the metallicity vary with overden- 
sity? Does the metallicity vary with redshift? How much scat- 
ter is there in the metallicity for a fixed density and redshift? 
Does this scatter change with density or redshift? Are the 
metals photoionized or does collisional ionization dominate? 
How does the metallicity change with the assumed spectral 
shape of the UV-background? What are the relative abun- 
dances of the different elements? Except for the last, which 
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we will address in future publications, this paper will address 
all of these questions. 

We measure the distribution of carbon in the IGM by apply- 
ing an extension of the pixel optical depth method, which we 
developed and tested in Aguirre, Schaye, & Theuns (2002, 
hereafter Paper I) — which was itself an extension of the 
earlier work by Cowie & Songaila (1998) (see also Songaila 
1998, Ellison et al. 2000, and Schaye et al. 2000a) — to a set 
of high-quality quasar spectra taken with the Keck and VLT 
telescopes. The basic idea is to measure the distribution of 
C IV pixel optical depths as a function of the corresponding 
H I optical depth, corrected for noise and absorption by other 
transitions, and to convert this correlation into an estimate of 
the metallicity as a function of the density. Our method differs 
in two respects from that proposed in Paper I: (1) we measure 
the full distribution of metals instead of just the median metal- 
licity and (2) we correct for noise and contamination using the 
data themselves instead of simulations. Furthermore, we use 
the C III/C IV ratio to set an upper limit on the temperature of 
the enriched gas, ruling out collisional ionization as the dom- 
inant ionization mechanism. 

This paper is organized as follows. In ^2]we describe our 
sample of quasar spectra, and in ^3]we use this sample to mea- 
sure the evolution of the mean absorption, which we need to 
normalize the intensity of our U V b ackground models. Be- 
fore describing these models i n 34.21 we briefly describe our 
hydrodynamical simulation in 34.11 The simulation is used to 
compute the gas density and temperature as a function of the 
H I optical depth, which are needed to com pute the ionization 
correction factors, as discussed in detail in 34.31 In 5j5]we pro- 
vide a step-by-step descr iptio n of our method for measuring 
the median metalli city ( 35. 1> and the distribution of metals 
at a fixed density ( 35.2> . To make the reading more interest- 
ing, we illustrate the method by presenting results for one of 
our best quasar spectra. Section[6|contains our measurements 
of the C III/C IV ratio, which support our assumption that the 
temperature is low enough for collisional ionization of C IV 
to be unimportant. In 0we present our measurements of the 
distributi on o f carbon as a function of overdensity and red- 
shift. In 37.1l we show how the results change if we vary the 
UV background. In ^8]we discuss the results, compute mean 
metallicities and filling factors, compare with previous work, 
and estimate the size of systematic errors. Finally, we sum- 
marize our conclusions in |9] 

How should this paper be read? Given the length of the pa- 
per, we encourage reading the conclusions (|9} first. Readers 
who are not interested in the details of the method would in 
addition only need to read §S0and[8] Those who would like 
to understand the method, but are not interested i n knowin g all 
of the details, would benefit from also reading § 34.2ll4~3l and 
13 Readers who would like to know more about the mean ab- 
sorption or the constraints on the ionization mechanism may 
want to read §|3]and|6] respectively. 

2. OBSERVATIONS 

We analyzed spectra of the 19 quasars listed in the first col- 
umn of Table [2 Fourteen spectra were taken with the UV- 
Visual Echelle Spectrograph (UVES; D'Odorico et al. 2000) 
on the VLT and five were taken with the High Resolution 
Echelle Spectrograph (HIRES; Vogt et al. 1994) on the Keck 
telescope (see col. [7] of Table [Q. The UVES spectra were 
taken from the ESO archive. Only spectra that have a signal- 
to-noise ratio (S/N) > 40 in the Lya region and that were 
publicly available as of 2003 January 31 were used. The 



UVES spectra were reduced with the ESO-maintained MI- 
DAS ECHELLE package (see Kim, Cristiani, & D'Odorico 
2001 and Kim et al. 2003 for details on the data reduction). 
The reduction procedures for the HIRES spectra are described 
in Barlow & Sargent (1997). The spectra have a nominal ve- 
locity resolution of 6.6 kms" 1 (FWHM) and a pixel size of 
0.04 and 0.05 A for the HIRES and UVES data, respectively. 

For a quasar at redshift z em (the emission red- 
shift), we analyze data in the redshift range 

(1 + Zem)A L y,3/ALy Q ~ 1 < Z < Zem ~ (1 + Zem)Av/c, where 

Av = max(4000,8Mpci/(z)//z) kms" 1 and H(z) is the Hubble 
parameter at redshift z extrapolated from its present value 
(H = 100ft kms-'Mpc" 1 ) assuming (f2 m ,fi A ) = (0.3,0.7). 
The lower limit ensures that Lya falls redwards of the 
quasar's Ly/3 emission line, thereby avoiding confusion with 
the Ly/3 forest. The region close to the quasar is excluded to 
avoid proximity effects. The minimum, median, and maxi- 
mum absorption redshifts considered are listed in columns 
(3), (4), and (5) of Table[T] 

Regions thought to be contaminated by absorption features 
that are not present in our simulated spectra were excluded. 
Examples are absorption by metal lines other than C III, C IV, 
Si III, Si IV, O VI, N V, and Fe II, and atmospheric lines. Seven 
spectra (Ql 101-264, Ql 107+485, Q0420-388, Q1425+604, 
Q2126-158, Q1055+461, and Q2237-061) contain Lya lines 
with damping wings (i.e., V(HI) > 10 19 cm" 2 ). Regions con- 
taminated by these lines or their corresponding higher or- 
der Lyman series lines were also excluded. For the case of 
damped Lya lines the excluded regions can be large (up to 
several tens of angstroms). Contaminating metal lines were 
identified by eye by searching for absorption features corre- 
sponding to the redshifts of strong H I lines and to the red- 
shifts of the absorption systems redward of the quasar's Lya 
emission line (trying a large number of possible identifica- 
tions). Due to the robustness of our method, which uses the 
C IV doublet to correct for contamination and is based on non- 
parametric statistics, our results are nearly identical if we do 
not remove any of the contaminating metal lines. 

3. THE OBSERVED EVOLUTION OF THE MEAN LYQ ABSORPTION 

Our measurements of the effective optical depth, r e ir = 
(— InF), where F is the normalized flux in the Lya forest, 
are pl otte d as a function of redshift in Fig.^and are listed in 
Table lB 1 l in Appendix B. Besides the quasars listed in Tabled 
we measured r e ff for two additional HIRES quasar spectra 
from the Rauch et al. (1997) sample: Q1442+293 (z em = 2.67) 
and Q0000-262 ( em = 4.1 1). These spectra are not part of our 
main sample because they do not have sufficient coverage to 
be useful for our metal analysis. Each Lya forest spectrum 
was divided in two, giving 42 data points in total. The left- 
hand panel shows the results obtained using all pixels in the 
Lya forest range; the right-hand panel shows the results ob- 
tained after removing all pixels that are thought to be con- 
taminated by metal lines, H I lines with damping wings, and 
atmospheric lines. Comparison of the two panels shows that 
removing the contamination generally reduces both the mean 
absorption and the scatter in the absorption at a given redshift, 
particularly at lower redshifts. 

The dotted lines show least-squares power-law fits to the 
data points. Note that the scatter is much larger than expected 
given the error bars: the x 2 /dof is about 5.2 and 2.4 for the 
left- and right-hand panels respectively (for 40 dof). This 
probably indicates that most of the remaining scatter is due 
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Table 1. Observed quasars 
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to cosmic variance (see also Kim et al. 2002). 

The dashed curve in Figure \l\ indicates the evolution of 
the Lya optical depth measured by Bernardi et al. (2003), 
who applied a novel x 2 minimization technique to 1061 low- 
resolution QSO spectra (their S/N > 4 sample) drawn from 
the Sloan Digital Sky Survey database. The Bernardi et al. fit 
gives absorption systematically higher by about 0. 1 dex than 
we obtain from our high-resolution spectra. For z > 3 our lo- 
cal continuum fits may underestimate the absorption and the 
Bernardi et al. results may be more reliable. However, fol- 
lower redshifts the continuum is well-defined and our results 
may be more accurate than those of Bernardi et al., whose 
spectra have insufficient resolution to allow for the identifica- 
tion and removal of metal lines. 

4. SIMULATION AND IONIZATION CORRECTIONS 

Because the C IV optical depth can tell us only about the 
density of C IV ions, we need to know the fraction of carbon 
that is triply ionized to determine the density of carbon. The 
ionization balance depends on the gas density, temperature, 
and the ambient UV radiation field. We use (variants of) the 
models of Haardt & Madau (2001) for the UV background 
and a hydrodynamical simulation to compute the density and 
temperature as a function of the H I optical depth. Before 
explaining ho w we convert the C IV/H I optical depth into a 
metallicity in 34.31 we will briefly describe our si mula tion and 
our method for generating synthetic spectra in 34.11 and our 
models for the metagalactic UV radiation in 34.21 

4. 1 . Generation of synthetic spectra 

We use synthetic absorption spectra generated from a hy- 
drodynamical simulation for two purposes: (1) to determine 
the gas density and temperature (which are needed to com- 
pute the ionization corrections) as a function of the H I opti- 
cal depth and redshift and (2) to verify that simulations using 
the carbon distribution measured from the observed absorp- 
tion spectra do indeed reproduce the observed optical depth 
statistics. 

The simulation that we use to generate spectra is identi- 
cal to the one used in Paper I, and we refer the reader to 



§3 of that paper for details. Briefly, we use a smoothed 
particle hydrodynamics code to model the evolution of 
a periodic, cubic region of a (n m ,n\,flt,h 2 ,h,a^,n,Y) = 
(0.3,0.7,0.019,0.65,0.9, 1.0, 0.24) universe ofcomoving size 
12/i _1 Mpc to redshift z = 1.5 using 256 3 particles for both the 
cold dark matter and the baryonic components. The gas is 
photoionized and photoheated by a model of the UV back- 
ground, designed to match the temperature-density relation 
measured by Schaye et al. (2000b). Since we recalculate 
the ionization balance of the gas when computing absorption 
spectra, this choice of UV background only affects the ther- 
mal state of the gas. 

The software used to generate the simulated spectra is also 
described in detail in Paper I (§3). Briefly, many sight lines 
through different snapshots of the simulation box are patched 
together 9 to create long sight lines spanning z = 1 .5 to z em - Ab- 
sorption from Lyl (A1216), Ly2 (A1026), . . . , Ly31 (A913), 
C III (A977), C IV (AA1548,1551), N V (AA1239,1243), O VI 
(AA1032,1038), Si III (A1207), Si IV (AA1394,1403), and Fe II 
(A1145, A1608, A1063, A1097, A1261, A1122, A1082, A1143, 
Al 125) is included. The long spectra are then processed to 
match the characteristics of the observed spectrum they are 
compared with. First, they are convolved with a Gaussian 
with FWHM of 6.6 kms -1 to mimic instrumental broadening. 
Second, they are resampled onto pixels of the same size as 
were used in the observations. Third, noise is added to each 
pixel. The noise is assumed to be Gaussian with a variance 
that has the same dependence on wavelength and flux as the 
noise in the observations. 

The ionization balance of each gas particle is computed 
from interpolation tables generated using the publicly avail- 
able photoionization package CLOUDY 10 (ver. 94; see Fer- 
land et al. 1998 and Ferland 2000 for details), assuming the 
gas to be optically thin and using specific models for the UV 
background radiation, which we will discuss next. 

9 Unlike in Paper I, we do not cycle the short sight lines periodically as this 
could potentially create spurious correlations between transitions separated 
by approximately an integral number of box sizes. 

10 See http : / /www . pa . uky . edu/~gary /cloudy. 
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FIG. 1 . — Observed effective optical depth as a function of redshift for all pixels (left-hand panel) and after removal of pixels contaminated by metal lines, Lya 
lines with damping wings, or atmospheric lines (right-hand panel). Horizontal error bars indicate redshift ranges used, vertical error bars are lcr errors measured 
by bootstrap resampling the spectra using chunks of 100 pixels. The solid curve shows the average effective optical depth in the simulated spectra; the dotted 
lines are power-law least-squares fits to the data points: logr c ff = log to + alog[(l + z)/4]. The best-fit values are (log to, a) = (—0.40 ± 0.02,3.40 ± 0.23) and 
(-0.44 ± 0.01,3.57 ± 0.20) for the left- and right-hand panels, respectively. The errors on the parameters of the fit were estimated by bootstrap resampling the 
data points. The dashed curve indicates the evolution of the Lya optical depth measured by Bernardi et al. (2003) using a \ 2 minimization technique on QSO 
spectra (their S/N > 4 sample) drawn from the Sloan Digital Sky Survey database. 



4.2. UV background models 

We use the models of Haardt & Madau (2001, hereafter 
HMO 1 ) 1 1 for the spectral shape of the meta-galactic UV/X-ray 
background radiation. Our fiducial model, which we will re- 
fer to as model QG, includes contributions from both quasars 
and galaxies, but we will also compute some results for model 
Q which includes quasars only (this model is an updated ver- 
sion of the widely used Haardt & Madau 1996 model). Both 
models take reprocessing by the IGM into account. 

The spectra have breaks at 1 and 4 ryd as well as peaks due 
to Lya emission. For model QG the spectral index -dJ^/dv, 
where J v is the specific intensity at frequency v, is about 1 .5 
from 1 to 4 ryd, but above 4 ryd the spectrum hardens consid- 
erably to a spectral index ranging from about 0.8 at z ~ 2 to 
about 0.3 at z ~ 4. For model QG (Q) the softness parameter 
S = Tei/Teen where T,- is the ionization rate for element i, in- 
creases from 350 (175) at z ~ 2, to 600 (230) at z ~ 3, to 950 
(280) at z w 4. 

To see how the results would change if the UV background 
were much softer at high redshift, as may be appropriate dur- 
ing and before the reionization of He II, we also create model 
QGS, which is identical to model QG except that the intensity 
above 4 ryd has been reduced by a factor of 10. 

All UV background spectra are normalized (i.e., the inten- 
sities are multiplied by redshift-dependent factors) 12 so that 
the effective optical depth in the simulated absorption spec- 
tra (Fig. [2 solid curves) reproduces the evolution of the ob- 
served mean transmission. We find that our simulation agrees 
with the observations if the H I ionization rate Thi/IO -13 pa 
8.7, 5.4, and 3.6 at z = 2, 3, and 4, respectively (note that 
since V oc O^/i 3 [e.g., Rauch et al. 1997], these values would 
be about 27% higher for the currently favored cosmology 
[n b h 2 = 0.0224, h = 0.71; Spergel et al. 2003]). 

1 1 The data and a description of the input parameters can be found at 

http://pitto. mib . inf n . it /^haardt /ref mo del .html. 

12 Multiplying the QG spectra by 0.58 for z < 2.5 and 0.47 for z > 2.5 
gives satisfactory results. 



4.3. Ionization corrections 

The ratio tqv /thi is proportional to the ratio of the C IV 
and H I number densities of the gas that is responsible for 
the absorption at that redshift. Converting this into a carbon 
abundance requires an estimate of the fraction of carbon that 
is triply ionized and the fraction of hydrogen that is neutral. 
These fractions depend on the gas density, temperature, and 
the spectrum of the UV background radiation. 

Both (semi-)analytic models (e.g., Bi & Davidsen 1997; 
Schaye 2001) and hydrodynamical simulations (e.g., Croft et 
al. 1997; Zhang et al. 1998; Paper I) suggest that there is a 
tight correlation between the H I Lya optical depth (or column 
density) and the gas density. For a given density thi depends 
weakly on the temperature and is inversely proportional to the 
H I ionization rate (assuming ionization equilibrium). 

We use our hydrodynamical simulation, which predicts a 
thermal evolution consistent with the measurements of the 
IGM temperature of Schaye et al. (2000b), to compute inter- 
polation tables of the density and temperature as a fun ction 
of the Lya optical depth and redshift. As discussed in 34.21 
we rescale the UV background so that the simulated spectra 
reproduce the observed mean Lya transmission. 

Because the absorption takes place in redshift space, multi- 
ple gas elements can contribute to the optical depth of a given 
pixel. We define the density/temperature in each pixel as the 
average over all gas elements that contribute to the absorp- 
tion in the pixel, weighted by their H I optical depths (Schaye 
et al. 1999). Therefore, the densities we quote are effectively 
smoothed on the same scale as the Lya forest spectra. 

Thermal broadening, the differential Hubble flow across 
the absorbers, and peculiar velocity gradients all contribute 
to the line widths. Simulations show that the first two are 
typically most important for the low column density for- 
est (e.g., Theuns, Schaye, & Haehnelt 2000). The line 
widths (^-parameters) are typically ~ 20-30 kms -1 for 
H I (e.g., Carswell et al. 1984) and - 5 - 15 kms" 1 for 
(directly detectable) C IV (e.g., Rauch et al. 1996; The- 
uns et al. 2002b). The thermal broadening width is b± = 
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13 kms" 1 (m H /m) I;/2 (7y 10 4 K) 1 / 2 , where m is the atomic 
weight (me ph 12rae). The real-space smoothing scale is sim- 
ilar to the local Jeans scale, which varies from ~ 10 2 kpc 
for the low column density forest to ~ 10 kpc for the rare, 
strong absorption lines arising in collapsed halos (e.g., Schaye 
2001). For our cosmology, a velocity difference Av corre- 
sponds to a physical scale of 34 kpc (Av/10 kms" 1 ) at z = 3. 
Thus, our densities are typically smoothed on a scale of 50 - 
100 kpc. 

When Hubble broadening dominates, as is the case for all 
but the highest overdensities studied in this work, there is little 
ambiguity in the relation between optical depth and gas den- 
sity. However, if thermal broadening were to dominate, then 
the absorption in pixels with low optical depths could in prin- 
ciple arise in the thermal tails of high-density gas. This would 
introduce scatter in the relation between H I optical depth and 
overdensity. Furthermore, because carbon is heavier than hy- 
drogen, it would result in discrepancies between the densities 
of the gas responsible for the C IV and H I absorption at a fixed 
redshift. 13 

Fortunately, simulations indicate that there is little scatter 
in the density corresponding to a fixed optical depth. Fig. 6 of 
Paper I shows that the relation between thi and 8 remains tight 
out to the maximum measurable optical depths, thi ~ 10 2 , 
while Fig. 7 of Paper I shows that weighting the gas by rciv 
and thi gives nearly the same densities, although the differ- 
ence does increase with overdensity. 

We can understand the empirical result that the relation be- 
tween gas density and optical depth remains well-defined up 
to (at least) S ~ 10 2 as follows. First, absorbers corresponding 
to overdensities 6 < 10 2 (Nm < 10 16 5 cm" 2 at z = 3) are gen- 
erally not isolated, compact clouds. The absorption arises in a 
fairly smooth gas distribution, which often contains multiple 
peaks. The substructure is generally resolved in redshift space 
because the differential Hubble broadening across the absorp- 
tion complex and/or internal peculiar velocity gradients are 
not much smaller than the thermal broadening scale (which is 
in turn greater than the instrumental resolution for H I). More 
importantly, because the column density distributions of H I 
and C IV are steep, most low optical depth pixels are not near 
high optical depth pixels; i.e., the contribution of "local" gas 
to the optical depth in any given pixel is generally greater than 
that of the thermal tails of high column density absorbers. 

In Paper I we used H I weighted quantities to compute the 
H I fraction and C IV weighted quantities to compute the C IV 
fraction. Since making this distinction changes the results 
only marginally (and only for the highest densities, see Fig. 7 
of Paper I) and since it is not clear what density to assign the 
resulting metallicity to, we use only H I weighted quantities 
here. 

Figure|2]shows contour plots of the overdensity (left panel) 
and temperature (right-hand panel) as a function of thi and z- 
Note that the hydrogen number density corresponding to an 

13 The fact that the CIV and HI fractions scale differently with the gas 
density and temperature also contributes to the differences between the thi 
and tciv weighted densities. However, we find from our simulation that the 
difference in the atomic weights is more important, except perhaps if the CIV 
and HI absorption arises in different gas phases. Note that the effects of the 
difference in the thermal broadening scales of HI and CIV could in principle 
be partially compensated for by smoothing the spectra on a scale greater than 
^th.civ, but smaller than foth.Hl- However, we find that doing this improves the 
results only marginally for CIV. 



overdensity 5 is given by 

Figure [3] shows the ionization correction as a function of 
temperature and density for the UV-background models QG 
(solid contours) and Q (dashed contours), all for z = 3. The 
ionization correction is defined as the factor with which the 
optical depth ratio tciv / Tin must be multiplied to obtain the 
carbon abundance relative to solar (the contours are labeled 
with the log of this factor). Note that the ionization correction 
can be modest even when both the H I and C IV fractions are 
negligible (as is the case for T ^> 10 5 K). 

The shape of the contours in Fig.[5]is easily understood. At 
high temperatures, collisional ionization dominates and the 
ionization balance becomes independent of the density. At 
lower temperatures, the change in the ionization correction 
with density can be accounted for by the degree to which car- 
bon is ionized. At high densities, the correction is large be- 
cause most of the carbon is less than triply ionized, and at 
very low densities, it is large because most of it is more than 
triply ionized. Compared with model QG, model Q predicts 
a lower C IV fraction (and a higher ionization correction) for 
low-density gas because it has more photons that can ionize 
C IV (which has an ionization potential of about 4.7 ryd). The 
difference between the spectra of models QG and Q is par- 
ticularly large above 4 ryd because stars produce very few 
photons above the He II Lyman limit compared with quasars 
(recall that the models have been scaled so that they have iden- 
tical H I ionization rates). 

For both models the C IV fraction is near a maximum in the 
density range where our data are best, -5 < log «h ^$ _ 4, and 
consequently the ionization correction is relatively insensitive 
to the density over the range that is of interest to us. The 
correction is also insensitive to the temperature as long as the 
gas is predominantly photoionized. In fact, for -5 < log«H ;$ 
-4, the temperature dependence is weak even for temperatures 
as high as 10 6 K. 

5. METHOD AND RESULTS FOR Q1422+230 

In this section we will discuss our method for measuring the 
distribution of metals in the diffuse IGM. We will illustrate the 
method by showing results for one of our best quasar spectra: 
Q1422+230. 

5.1. The median metallicity as a function of density 

Our method for measuring median metal abundances as 
a function of density from QSO absorption spectra is pre- 
sented in detail in Paper I. Readers unfamiliar with pixel opti- 
cal depth techniques may benefit from reading §2 of Paper I, 
which contains a general overview of the method. 

Our goal is to measure the abundance of carbon as a func- 
tion of the gas density and redshift, but our observable is the 
quasar flux as a function of wavelength. In order to derive the 
metallicity we need to do the following: 

1. Normalize the spectra. The quasar spectra are divided by 
a continuum fit to obtain the normalized flux F = F \, s /F cont 
as a function of wavelength, where F Q \, S and F cont are the ob- 
served and continuum flux respectively. The noise array is 
similarly rescaled: u = cr b s /F cont . The HIRES and UVES 
spectra were continuum fitted as described in Rauch et al. 
(1997) and Kim et al. (2001), respectively. 

We found that relative to their signal-to-noise ratio most 
UVES spectra had larger continuum fitting errors than the 
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FIG. 2. — Left: Contour plot of log 8 as a function of log thi and redshift z. Overdensity increases with optical depth but decreases with redshift. Right: Contour 
plot of log T as a function of log thi and z. The temperature is relatively constant over the parameter range of interest. 
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log T (K) 

FIG. 3. — Ionization correction factor, [C/H] — logfrciv/THi), as a function 
of the temperature and the hydrogen number density. Solid (dashed) contours 
are for the UV background model QG (Q). At high temperatures collisional 
ionization dominates and the ionization balance becomes independent of the 
density. The temperature below which photoionization becomes important 
(and the contours start to curve) is higher if the density is lower. The results 
are only sensitive to the hardness of the UV-background if the density is low. 
In the regime of interest to us (n H ~ 10~ 5 to lfT 4 cm~ 3 , T ~ 10 4 - 10 4 5 K), 
the ionization correction is near a minimum and thus relatively insensitive to 
the exact density and temperature. 



HIRES spectra. All C IV regions were therefore renormal- 
ized using the automatic procedure described in §4. 1 of Paper 
I: We divide the spectra into bins of rest-frame size 20 A with 
centers and find for each bin the median flux fa. We then 
interpolate a spline across the fi and flag all pixels that are 
at least N z Ja below the interpolation. The medians are then 
recomputed using the unflagged pixels and the procedure is 
repeated until the fit converges. Finally, the flux and error ar- 
rays are rescaled using the fitted continuum. Tests using sim- 
ulated spectra indicate that N~ = 2 is close to optimal for the 
C IV region, giving errors in the continuum that are smaller 
than the noise by an order of magnitude or more. When ap- 
plied to the observations we find that the adjustments in the 
continuum are typically smaller than the noise by a factor of 
a few for the UVES spectra and by an order of magnitude for 



the HIRES spectra. 

2. Recover the H I and C IV optical depths as a function of 
redshift. The recovery r = -InF is imperfect because of con- 
tamination by other absorption lines, noise, and continuum 
fitting errors. In Paper I (§4 and Appendix A) we describe 
methods 14 to partially correct for these effects, such as the use 
of higher order Lyman series lines to estimate the H I optical 
depth of saturated absorption features [i.e., F(X) < 3cr(A)] and 
an iterative procedure involving both components of the C IV 
doublet to correct for (self-)contamination. We will refer to 
the values resulting from these procedures as "recovered" op- 
tical depths. 

Pixels for which the H I optical depth cannot be accurately 
determined because all available Lyman series lines are satu- 
rated are removed from the sample. Pixels that have negative 
C IV optical depth (which can happen as a result of noise and 
continuum fitting errors) are given a very low positive value 
in order to keep the ranking of optical depths (and thus the 
medians) unchanged. 

3. Bin the pixels according to thi and compute the me- 
dian rciv for each bin. There are two reasons for binning in 
thi. First, binning allows us to use nonparametric statistics, 
such as the median rciv and other percentiles, which makes 
the method much more robust. This is important because our 
optical depth recovery is imperfect. Second, the H I optical 
dept h is t hought to be correlated with the local gas density 
(see 34. 3> . which is needed to compute the ionization balance. 

The errors in the medians are computed as follows. First, 
the Lya forest region of the spectrum is divided into N chunks 
of 5 A, and H I bins which contain less than 25 pixels or con- 
tributions from fewer than five chunks are discarded. Sec- 
ond, a new realization of the sample of pixel redshifts is con- 
structed by bootstrap resampling the spectra; i.e., N chunks 
are picked at random with replacement and the median tciv 
is computed for each H I bin. This procedure is repeated 100 
times. Finally, for each H I bin we take the standard devia- 
tions of the median log rciv as our best estimate of the errors 

14 Our method differs from that of Paper I in our treatment of saturated 
metal pixels, i.e., CIV pixels with a normalized flux F(A) < 3<r(A), where 
u(A) is the normalized noise array. Instead of setting the optical depth in 
these pixels to — ln3a"(A) as was done in Paper I, we set it to a much larger 
value (10 4 ). 
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in the medians (the distribution of the medians from the boot- 
strap realizations is approximately lognormal, although with 
somewhat more extended tails). 

The top left-hand panel of Figure |4] shows the results for 
Q1422+230, one of our best quasar spectra. The recovered 
C IV and H I optical depths are correlated down to logTm « 
0.1, below which logrciv ~ -3.0 independent of thi. In Paper 
I we demonstrated that the correlation flattens at low optical 
depth due to the combined effects of noise, contamination, 
and continuum fitting errors. 

A similar analysis was applied to an even higher S/N spec- 
trum of the same quasar by Ellison et al. (2000). Our results 
agree with theirs for log thi > 0.3, below which Ellison et al.'s 
correlation flattens off at logrciv ~ -2.8. We thus detect the 
Tcrv - t hi correlation down to smaller optical depths despite 
the fact that our data has an S/N about a factor of 2 lower 
than the spectrum of Ellison et al. (2000). The improvement 
is mainly due to our correction for (self-)contamination. 

Note that there is additional information contained in the 
distribution of tciv that is not used when we only consider 
the median. Different from Paper I, we will also consider 
other percentiles than the 50th (i.e., the median), which will 
allow us to characterize the scatter in the metallicity at a fixed 
densi ty. Our method to measure this scatter is discussed in 

3221 

4. Correct the median tciv(thi). The correlation between 
the recovered tcw and thi flattens off at low thi where noise, 
continuum fitting errors, and contamination wash out the sig- 
nal. In Paper I we used simulations to calibrate the difference 
between the true and recovered tcw and then used this in- 
formation to correct the tctv- Because we would like to mini- 
mize our use of the simulations and because the correction for 
C IV does not show any dependence on thi, we choose to cor- 
rect the recovered C IV optical depth using the observations 
themselves. 

Assuming that, as is the case in the simulations (see Pa- 
per I), the asymptotic flat level tciv(thi — > — oo) is the median 
signal due to contamination, noise, and/or continuum fitting 
errors and that this spurious signal is independent of thi, we 
can correct for this component to the signal by subtracting it 
from the data points. 

The asymptotic flat level is determined as follows. We first 
find r c , the H I optical depth below which the correlation van- 
ishes, by fitting a power law to the data 

log(rciv)={^ + Blog(THi/Tc) 2tr c (2) 

Next, we define the asymptotic C IV optical depth, r m i n , as the 
median tciv of those pixels that have thi < r c , 

median(rciv|T H i < r c ) (3) 

(and similarly for other percentiles than the median). The hor- 
izontal and vertical dotted lines in the top left-hand panel of 
Figure|4]indicate r m i n and r c , respectively. As was the case for 
the data points, the error on t^,, was computed by bootstrap 
resampling the spectrum. To be conservative, we add the error 
in Tmin linearly (as opposed to in quadrature) to the errors on 
the recovered tciv- Finally, for display purposes only, we fit 
equation Q again to the data, but this time keeping r m i„ fixed 
at the value obtained before. The result is shown as the solid 
curve in the figure. 

The corrected data points are shown in the top right-hand 
panel of Figure |4] and should be compared with the dashed 
curve in the same panel, which connects the recovered data 



points that were plotted in the top left-hand panel. The cor- 
rection is significant only for those points that are close to 
Tmin- Data points with thi < r e are effectively converted into 
upper limits. 

5. Compute the carbon abundance as a function of density. 

Assuming ionization equilibrium, which should be a very 
good approximation for gas photoionized by the UV back- 
ground, we can compute the ionization balance of hydrogen 
and carbon given the density, temperature, and a mode l for the 
ionizing background radiation. As discussed in 34. 31 we use 
CLOUDY to compute ionization balance interpolation tables 
as a function of redshift, density, and temperature, and we use 
our hydrodynamical simulation to create tables of the density 
and temperature as a function of redshift and H I optical depth. 

To convert the optical depth measurements into metallic- 
ity estimates, we first compute the median redshift and thi 
of each H I bin and then use the interpolation tables to com- 
pute the corresponding density, temperature, and ionization 
balance. For example, the dashed curve in the bottom left 
panel of Figure|4]shows C/C IV, the inverse of the C IV frac- 
tion, as a function of thi. The solid curve shows the ionization 
correction factor: 10 [C / H] /(tciv/7"hi)- 

Given the corrected C IV optical depth rcrv> the recovered 
H I optical depth thi, and the ionization correction, we can 
compute the carbon metallicity: 

[C/HI ^(S^^™)-' C/H '- <4 > 

where /; and A, are the oscillator strength and rest wavelength 
of transition i, respectively (/crv = 0.1908, f m = 0.4164, 
Aciv = 1548.2041 A, A ffl = 1215.6701 A), and we use the so- 
lar abundance (C/H)q =-3.45 (number density relative to hy- 
drogen; Anders & Grevesse 1989). The resulting data points 
are plotted in the bottom right-hand panel of Figure |^ which 
shows the carbon metallicity as a function of the overdensity 
(bottom axis) or Lya optical depth (top axis). 

Metals are detected over two decades in density. The me- 
dian metallicity increases from a few times 10~ 4 around the 
mean density to ~ 10~ 3 to 10~ 2 around an overdensity S ~ 10 2 . 
The best-fit 15 constant metallicity is [C/H] = -3.30^;g|, but 
this fit is ruled out at high confidence (% 2 /dof ps 30.6/9 w 3.4 
for the data points with r H i > t c , corresponding to a prob- 
ability Q = 3.5 x 10 -4 ). The best-fit power-law metallic- 
ity is [C/H] =-3.12^ + 0.90^;J|(log5-1.0), which has 
X 2 /dof « 3.42/8 re 0.43, or a probability Q = 0.91. 

The errors in this plot reflect only the errors in the cor- 
rected median rciv (top right-hand panel) and do not take 
uncertainties in the ionization balance into account. The dot- 
ted (dashed) curves illustrate how the metallicity changes if 
the temperature (density) is changed by a factor of 2. The 
metallicity is most sensitive to the assumed density for high 
5, where the ionization correction factor changes rapidly with 
thi. Fortunately, the ionization correction is relatively insen- 
sitive to both the density and the temperature in the regime 
where our data is best. 

5 All fits in this paper are obtained by minimizing x 2 - 
(K '— Z(5j)]/ci) > where Z, = [C/H], are the metallicity data points and 
Z(<5) is the function being fitted. The sum is over all data points for which 
thi > To The errors in the recovered tciv are lognormal, but the subtraction 
of T m j„ results in asymmetric error bars for [C/H]. We therefore compute 
a closely spaced grid of errors (e.g., ±0.01ct,0.02(t, . . .) in tciv and T ra j„, 
propagate these to obtain an equivalent (irregular) grid of errors for [C/H], 
and then compute \ 2 using the latter grid. 
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FIG. 4. — Measuring the carbon abundance in Q1422+230. Data points are plotted with 1 and 2 cr error bars. Top left: Recovered tciv1548 plotted as a function 
of the recovered 7hh216- For thi < r c (vertical dotted line) the signal is lost due to contamination, noise, and/or continuum errors and tciv ~ T m j„ (horizontal 
dotted line). Top right: Corrected tciv vs. recovered thi. Data points have been corrected by subtracting T m j„ from the recovered tciv plotted in the top left-hand 
panel and indicated here by the dashed curve. The correction affects points close to T m i n , converting the lowest points into upper limits. Bottom left: Solid and 
dashed curves show the ionization correction factor and the inverse of the carbon fraction, respectively. The ionization correction is insensitive to the Lya optical 
depth in the regime where tciv is best measured. Bottom right: The carbon metallicity as a function of the overdensity (bottom axis) or thi (top axis). Data 
points were obtained by inserting the corrected tciv, me recovered thi (i.e., the data points shown in the top right-hand panel), and the ionization correction 
factor (solid curve in the bottom left-hand panel) into eq.j4j. The solid line shows the best-fit power law: [C/H] = — +O.90^j'[|(log<5 — 1.0). Dashed and 
dotted curves indicate the result of changing the density and temperature, respectively, by a factor of 2. 



5.2. The distribution of metals at a fixed density 

Having measured the median carbon abundance as a func- 
tion of the density, we use the simulation to check whether 
our measured metallicity profile can reproduce the observed 
optical depth statistics. The bottom set of data points in the 
left-panel of Figure [5] shows the observed median C IV opti- 
cal depth as a function of thi. The bottom solid curve shows 
the median tciv in the simulation (averaged over 10 spectra) 
for the best fit metallicity profile measured from the obser- 
vations [each gas particle was given a metallicity [C/H] = 
-3.12 + 0.90(log5- 1.0)]. Clearly, the simulation curve is a 
good fit of the observations (\ 2 probability 16 Q = 0.21), which 
gives us confidence that, given a model for the UV back- 

16 Because we do not want to compare simulated and observed noise, the 
X 2 is computed only for the data points with thi > T c , where r c is computed 
for the observed spectra. The light data points were excluded from the x 2 



ground, our method for measuring the median carbon abun- 
dance as a function of density is robust. 

Although the median tciv is reproduced by the simulation, 
this is not the case for other percentiles. From top to bottom 
the sets of data points in Figure|3]correspond to the 84th, 69th, 
and 50th percentiles. The simulation clearly cannot reproduce 
the width of the observed metallicity distribution (Q < 10~ 4 
for the higher percentiles). This indicates that the scatter in 
the observations is greater than can be accounted for by con- 
tamination and noise effects, both of which are present in the 
simulations. The scatter is also unlikely to be due to con- 
tinuum fitting errors because those would result in a scatter 
c(tciv) that is independent of thi, contrary to what is ob- 
served. Hence, to reproduce the observed optical depth statis- 

calculation. For the same reason, we have added T m i n o b s — T m i n s i m to the sim- 
ulated data points. The horizontal dotted lines indicate the original T m ; n s i m . 
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Fig. 5. — Comparison of the optical depth statistics of observed and simulated spectra using the metal distribution measured from the observations. From top 
to bottom the three sets of data points are the 84th (open diamonds), 69th (solid diamonds) and 50th (triangles) percentiles of the recovered CIV optical depth as 
a function of rm for Q1422+230. For clarity, the 84th and 69th percentiles have been offset by +1.0 and +0.5 dex, respectively. The curves in the left-hand panel 
are for a simulation in which each particle was given the median metallicity measured from the observations, [C/H] = -3.12 + 0.90(log<5- 1.0). The simulation 
can fit the observed median tqiv (x 2 probability Q = 0.21), but not the observed tciv(thi) f° r the other percentiles (Q < ICr^). The curves in the right-hand 
panel are for a simulation that has the same median metallicity, but which includes scatter. The simulation cube was divided into 10 3 cubic sections, and all 
particles in each section were given a metallicity of [C/H] = — 3.12+i+0.90(log 5— 1.0), where which is the same for all particles in the subvolume, is drawn 
at random from a lognormal distribution with mean and variance a = 0.81 dex as measured from the observations. The simulation provides an acceptable fit to 
all percentiles (from top to bottom Q = 0.33, 0.69, and 0.90). 



tics, we need to extend the method of Paper I to allow for 
scatter in the metallicity. 

Our method to reconstruct the median metallicity can of 
course also be applied to other percentiles. As we will show 
below, for a fixed density the metallicity distribution is ap- 
proximately lognormal. A lognormal distribution is parame- 
terized by the mean (logZ) = median(logZ) and the variance 
fi 2 (logZ), and the fraction / of pixels that have a value smaller 
than xa, is given by the Gaussian integral 

/W=-L/ e-' 2 / 2 dt. (5) 

For a given H I bin, we measure a by computing logZ for a 
large number of percentiles (i.e., values of /) corresponding 
to a regular grid of x-values, and fitting the linear function 
logZ = ax to the data points (we use a spacing Ax = 0.25 and 
require a minimum of four data points). The procedure is il- 
lustrated in Figure [6] which shows the measured metallicity 
as a function of x for three different densities. The data ap- 
pear well fitted by lognormal distributions {solid curves) with 
a w 0.97 for 6 « 0.9 (r ffl = 0.95; left-hand panel), a « 0.73 for 
5 w 4.3 (thi = 16.8; middle panel), and a « 0.64 for 5 s» 15.7 
(thi = 94.4, right-hand panel). Note that we cannot get reli- 
able errors on a from the fits shown in the figure because the 
errors on the data points are all correlated since the abscissa 
corresponds to a ranking. 
Figure [6] demonstrates that the distribution of pixel metal - 



licities is well described by a lognormal function over at least 
the range -0.5er to +2cr. To probe the distribution further into 
the low-metallicity tail we would need to reduce the noise, 
since the lowest percentile/x-value corresponds to the mini- 
mum positive Tcrv (the rest of the pixels have a flux greater 
than e" Tmi °). To probe further into the high-metallicity tail we 
would need more pixels since for pixels one cannot measure 
a percentile greater than (N— 1)/N [in practice we only make 
use of percentiles for which N(l— f) > 5 in order to prevent a 
few anomalous pixels from affecting the results]. 

The right-hand panel of Figure [5] compares the observed 
percentiles with a simulation that uses a lognormal metallic- 
ity distribution, with the same median metallicity as a function 
of density as before and with a constant scatter of cr([C/H]) = 
0.81 dex (the median value of the scatter measured in the ob- 
servations for the different densities) on a comoving scale 17 
of 1.2 h~ l Mpc. Contrary to the simulation without scatter 
{left-hand panel), this simulation is consistent with all the ob- 
served percentiles {\ 2 probability Q = 0.33, 0.69, and 0.90). 
To estimate the error in tr([C/H]) we compared the observed 
percentiles with those computed from simulated spectra us- 
ing varying amounts of scatter in the metallicity. We find that 
values of <r([C/H]) = 0.81 ±0.25 are acceptable. 

Unlike the median metallicity, the scatter is independent of 

17 The exact scale is unimportant as long as it is at least as large as the 
effective smoothing scale of the absorbers and smaller than the simulation 
box. 
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FIG. 6. — Distribution of carbon metallicities for three gas overdensities: S 1.6 (left), <5 sa 5.4 (middle), and 8 « 19.4 (right). For each density (i.e., HI bin) 
the CIV pixels were ranked accor ding to tciv and a number of percentiles were computed. The metallicity corresponding to each percentile was computed as 
described for the median in il5.ll The missing percentiles correspond to negative optical depths (i.e., the normalized flux is greater than unity). The curves are 
lognormal fits to the data points. The best-fit slopes are ct([C/H]) = 0.97, 0.73, and 0.64 for the left, middle, and right panels, respectively. 



the assumed spectral shape of the UV background radiation. 
This is because the ionization correction consists of multi- 
plying the observed optical depth ratio tciv /thi with a fac- 
tor that is independent of tcw'- [C/H] = log(rciv/THi)+^, 
where X is the log of the ionization correction factor (see eq. 
|4)) which depends only on thi- Hence, if log(Tcrv/7Hi) is 
distributed lognormally with mean /i and variance a 2 , then 
[C/H] is also distributed lognormally with variance a 2 , but 
with mean fj,+X. Note, however, that we have assumed the 
UV background to be uniform. If there are large fluctuations 
in the UV radiation at a fixed H I optical depth, then the true 
scatter in the metallicity could be smaller. 

Having found that the metallicity distribution is well fitted 
by a lognormal function, we can estimate the mean metallicity 
as follows 

log(Z)=log (—L^ [ +0 ° Ze< l ° sZ - {losZ) y /2a2 dlogz)(6) 
\v27rcr J-oo J 

= (logZ) + ^ ( i 2 «(logZ) + 1.156T 2 , (7) 

where (logZ) = median(logZ) and a are obtained as explained 
above. For Q 1422+230 this gives a mean metallicity of 
log (Z) re -2.36 ± 0.48 at 8 = 10. 

Unlike estimates based on the mean C IV optical depth, 
this estimate of the mean metallicity is based on a fit to the 
full metallicity distribution and is thus relatively insensitive 
to anomalous pixels. On the other hand, since we have only a 
finite number of pixels, this estimate of the mean does implic- 
itly assume that the metallicity distribution remains lognormal 
beyond where we can actually measure it. While the mean 
metallicity is insensitive to the shape of the low-metallicity 
tail of the distribution, the contribution of the unmeasured 
high-metallicity tail to the mean metallicity depends strongly 
on the variance. For example, cutting of the distribution at 
+2cr reduces the mean by factors of 2.6 and 1.2 for a = 1.0 
and 0.5 dex, respectively. For Q1422+230 we find a pa 0.8, 
which is small enough for the mean to be fairly insensitive to 
the unmeasured high-metallicity tail. 

5.3. Summary of results for Q1422+230 



Using model QG for the UV background and the 8(Tm,z) 
and T(tui,z) relations obtained from the simulations, we have 
found the following from analyzing the statistics of tciv(thi) 
inQ1422+230: 

1. The observed median optical depths are inconsistent 
with a constant metallicity. The observations are 
well fitted by a median carbon abundance that is a 
power law of the overdensity: [C/H] = -3.12^}'^ + 

0.90^;l|0og<S-1.0). 

2. The distribution of tciv at fixed thi is inconsistent with 
a metallicity that depends only on density. The obser- 
vations are well fitted by a lognormal metallicity distri- 
bution with a scatter of 0.81 ± 0.25 dex and a median 
that varies as above. 

6. THE CIII/CIV RATIO 

We compute the ionization corrections as a function of the 
H I optical depth and redshift using interpolation tables cre- 
ated from the simulation. Because the simulation reproduces 
the observed evolution of the mean absorption (Fig.^ and the 
temperature of the gas responsible for the low-column density 
Lya lines (Schaye et al. 2000b), we are confident that the ion- 
ization balance of the gas responsible for the H I absorption is 
on average well determined. However, it is possible in prin- 
ciple that the C IV and H I absorption at a given redshift is 
dominated by different gas parcels. For example, Theuns et 
al. (2002b) predict that many of the metals produced by dwarf 
galaxies reside in hot gas bubbles that do not contribute sig- 
nificantly to the H I absorption. In their simulation most of the 
carbon has T ~ 10 6 K and is too highly ionized to give rise to 
C IV absorption. In this scenario the observed C IV absorption 
arises in the low-temperature tail of the carbon temperature 
distribution. 

Although the ionization corrections are insensitive to small 
changes in the temperature (see Fig. [4}, they could be consid- 
erably in error if the temperature were high enough for colli- 
sional ionization to be important (T > 10 K for C IV). Fur- 
thermore, such a high temperature would change the relation 
between thi and 8, causing us to underestimate the density 
corresponding to a given H I bin. 
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FIG. 7. — Contour plot of log(Tcm/Tciv) in the density-temperature plane. 
Solid and dashed contours are for the z = 3 QG and Q UV backgrounds, 
respectively. For each density there is a critical temperature below which the 
CHI/CIV ratio is nearly constant and above which it decreases rapidly. For 
the densities of interest here («h ~ 10~ 5 to 10~ 3 ), the ratio varies from 0.1 to 
10. 



Fortunately, there is a way to constrain the temperature of 
the gas responsible for the carbon absorption by comparing 
the strengths of the C III (A977) and C IV (AA1548, 1551) 
transitions. Figure shows a contour plot of the predicted 
logTcm/rciv in the density-temperature plane. The solid 
(dashed) contours are for the z = 3 QG (Q) model of the UV 
background radiation. From the figure it can be seen that a 
measurement of the C III/C IV ratio yields an upper limit on 
the temperature independent of the density and the UV back- 
ground. 

The left-hand panel of Figure [8]shows the median C III op- 
tical depth as a function of tqw for Q 1422+230. The C IV 
optical depths were recovered as described in §5.11 while the 
C III optical depths were recovered as described in Paper I 
for O VI, i.e., contaminating higher order Lyman lines were 
removed. Absorption by C III is clearly detected over 1.5 
decades in tciv ■ For log tciv < - 1 .7 the apparent C III optical 
depth is constant, indicating that for these low C IV optical 
depths the true rem is smaller than the spurious signal result- 
ing from contamination by other absorption lines, noise, and 
continuum fitting errors. 

Note that log tciv = -1 .7 corresponds to a typical H I optical 
depth of a few tens (see Fig. 0J. Although there is consider- 
able scatter in the tciv corresponding to a fixed thi, this does 
suggest that we are detecting C III mostly in gas with den- 
sity <5 > 10. Indeed, the correlation between rem and thi (not 
plotted) is detectable only for log thi > 1 . 

We can see from Figure that if the carbon were purely 
photoionized (T -C 10 K), we would expect tciii/tciv val- 
ues between 1 and 10 at z ~ 3 and 6 > 10. The right-hand 
panel of Figure [8] shows the log of the ratio tchi/tciv as a 
function of rcrv> where the rem are the median values plotted 
in the left-hand panel but after subtraction of r m i n (horizontal 
dashed line). The flat level T m j n wa s determined in the same 
way as we did for tciv(thi) m §E3] (step 4). Percentiles other 
than the median (not plotted) give nearly identical Tcm/Vcrv 
ratios (after subtraction of the corresponding r m i„ percentiles), 
indicating that this ratio is rather uniform. For log tciv > -1 -5 
we detect C III and log rem /^civ ~ 0.5, as expected for pho- 



toionized gas at this redshift. Crucially, optical depth ratios 
typical for collisionally ionized C IV (T > 10 5 K) are ruled 
out by the data. Hence, the temperature of the gas responsible 
for the carbon absorption is smaller than 10 5 K, at least for 
densities <5 > 10. This provides an important constraint for 
theories of the enrichment of the IGM through galactic winds 
and gives us confidence that our assumption that photoioniza- 
tion dominates is reasonable. 

We have computed the C III/C IV ratio for all quasars in our 
sample that cover the C III region. The results are summa- 
rized in Figure|9] in which TcmAcrV; averaged over Tcrv> is 
plotted as a function of redshift. While this plot uses the me- 
dian rem, the other percentiles give nearly identical results. 
For all QSOs the ratio is in the range expected for photoion- 
ized gas with «h ~ 10~ 4 cm -3 and T < 10 5 K. The solid 
line shows the best linear least absolute deviation fit to the 
simulated versions of all the quasar spectra, using the metal- 
licity distributions measured from the observed spectra and 
the QG background. The dashed line shows a similar fit for 
the Q background. Overall, the simulations agree well with 
the observations, although they may slightly underpredict the 
C III/C IV ratios for z > 3. The QG background seems to fit 
the data slightly better than model Q, but the difference is too 
small to be significant. In the simulation the small decrease in 
the C III/C IV ratio with time is caused both by the decrease 
in the typical density of C IV absorbers due to the expansion 
of the universe and by the hardening of the UV background 
radiation. 

Thus, the observed C III/C IV ratios are in the range ex- 
pected for photoionized gas, but are inconsistent with scenar- 
ios in which a large fraction of the C IV absorption takes place 
in gas with temperatures T > 10 5 K. 

7. RESULTS FOR THE FULL SAMPLE 

In 5j5] we showed results for Q1422+230 to illustrate our 
method for measuring the distribution of carbon as a function 
of the gas density. In this section we will present the results 
from the complete sample of quasars. The results from the 
individual quasars are given in Appendix A. 

All abundances are by number relative to the total hydro- 
gen density, in units of the solar abundance [(C/H) Q = —3.45; 
Anders & Grevesse 1989]. 

Our goal is to measure the carbon abundance as a function 
of density and redshift. For a lognormal dist ribut ion, which 
we found to provide a good fit to the data (see 35. 2> . the distri- 
bution is determined by two parameters: the mean ([C/H]) = 
median([C/H]) and the standard deviation er([C/H]) (which 
we will often refer to as the scatter). Thus, we can character- 
ize the full distribution of carbon by fitting functions of the 
two variables 8 and z to all the data points for ([C/H] ) and 
er([C/H]) obtained from the individual quasars (see Fig. lAH . 
But before doing so, it is instructive to bin the data in one 
variable and plot it against the other. 

Figure Qol shows the median metallicity versus redshift for 
five different density bins of width 0.5 dex, centered about 
densities increasing from log (5 = -0.25 (top panel) to logS = 
1 .75 ( botto m right-hand panel). To facilitate comparison with 
Figure lATI we have plotted only one data point per quasar, ob- 
tained by fitting a constant metallicity to all the data points of 
the quasar that fall in the given density bin. The points have 
been labeled with the first four digits of the quasar name. The 
light data points have ler lower limits that extend to minus in- 
finity. The solid lines in each panel show the least -squ ares fits 
to the original data points (i.e., the ones from Fig. lAll not the 
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FIG. 8. — Tcm as a function of tqi\ for Q1422+230. Data points are plotted with 1 and 2a error bars. Left: Median recovered rem °f pixels binned as a 
function of recovered rciv- CHI is detected down to log rciv ~ —1.7. Right: The rcm/rciv ratio is plotted as a function of rciv- The median rem were corrected 
for noise, contamination, etc., as described for CIV in j|5. 1 K step 4). The ratio rcm/rciv is much higher than it would be if the temperature were higher than 
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FIG. 9. — The ratio rcm/rciv is plotted as a function of redshift for all 
Zem < 4 QSOs for which CHI is covered (we do not detect CIII for QSOs 
above this redshift). For each QSO the median rem were corrected for con- 
tamination, noise, etc., using the method described in J^lland the rcm/rciv 
ratio was averaged over all bins with rciv > T c (see 815.1 I for the definition 
of t c ). The solid line shows the best linear fit to the CIII/CIV ratio mea- 
sured from simulations of the same quasar spectra, which used the QG UV 
background and the metallicity distributions measured from the observations. 
The observed ratios agree well with the simulations, and comparison with 
Fig.l71shows that they are in the range expected for photoionized gas with 
n H ~ 1(T 4 ctrT 3 and T < 10 5 K. 



plotted rebinned data points) and the dotted curves indicate 
the la confidence limits. The errors on all fits were deter- 
mined by bootstrap resampling the quasars. Computing the 
errors using the \ 2 surface gives similar results, although the 
X 2 errors tend to be somewhat larger (smaller) than the boot- 
strap errors if the \ 2 P er degree of freedom is smaller (greater) 
than unity. The reduced \ 2 are somewhat small, particularly 
at low overdensities, indicating that we have overestimated 
the errors. This is not unexpected as we have used conserva- 
tive estimates of the errors on data points for which the correc- 
tion of th e "noise" component is significant (i.e., rciv ~ 7" m i n ; 
see EU- 



The first panel reveals that we have detected metals in un- 
derdense gas: [C/H] -3.6 at log 5 = -0.25, with a 2a lower 
limit of -4.12. We find that [C/H] > -5 at the 2Aa level 
(99.2% confidence). For percentiles higher than the 50th (i.e., 
the median) the detection is even more significant. For ex- 
ample, for the 69th percentile we find [C/H] > -5 at the 3.4cr 
level. Hence, barring errors in our estimate of the density con- 
trast, there is no question that a large fraction of underdense 
gas has been enriched. 

At the la level all bins are consistent with no evolution. 
The constraints are strongest for the bin log S = 0.5- 1.0, in 
which the maximum allowed increase in the median metallic- 
ity per unit redshift is 0.14 dex at the la level and 0.25 dex 
at the 2a level. Comparing the metallicity at z = 3 for the 
different overdensity bins, we see that the metallicity in- 
creases with overdensity from [C/H] ss -3.6 at log<5 -0.25 
to [C/H] w -2.6 at log i5w 1 .75. This trend may be more eas- 
ily analyzed using the next figure. 

Figure shows the median metallicity as a function of 
overdensity for three different redshift bins of width Az = 1.0 
centered on z = 2 , 3, and 4. The data points were taken directly 
from Figure IaTI As in the previous figure, light-colored data 
points have lcr lower limits of — oo. Note that many of the 
light-colored error bars correspond to data points that fall be- 
low the plotted range (some at — oo). Note that the data are 
not uniformly distributed in these redshift bins: the median 
redshifts of the data points in the three bins are z » 2.12, 2.84, 
and 4.07. Note also that the last bin contains only two quasars: 
Q1055+461 andQ2237-061. 

The three redshift bins give similar results, confirming that 
there is little evidence for evolution. We do, however, have a 
very significant detection of a positive gradient of metallicity 
with overdensity in the z = 2 a nd z = 3 bins (as seen also in 
the analysis of Q1422+230 in S15.H . The median metallicity 
increases from [C/H] « -4 at log 5 = -0.5 to [C/H] w -3 at 
log 5 = 1.0. The best-fit power law index a = d[C/H]/dlogS 
is 0.68^ ;}g for z = 2 and 0.72^ for z = 3. 

In addition to measuring a from the combined quasar data, 
we can also combine the a values measured in e ach quasar 
individually, which can all be read off from Figure lATl A lin- 
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FIG. 10. — Median metallicity as a function of redshift for different overdensity cuts. From the top left to the bottom right, panels correspond to the following 
density bins: log<5 = -0.5-0.0, 0.0-0.5, 0.5-1.0, 1.0-1.5, and 1.5-2.0. One data point per QSO is plotted, labeled with the first four digits of the QSO name. 
Horizontal error bars indicate redshift ranges, vertical error bars indicate 1 and 2<r errors. Light-colored data points have a lower l<r error of — oo. For each QSO 
and density bin the plotted data point was obtained by fitting a constant metallicity through all the HI bins of that QSO spectrum (plotted in Fig. lAll that fall in 
the density bin. This rebinning was done only for visualization purpo ses, facilitating direct comparison with Fig. lAll The least-squares fits, shown as the solid 
lines, were obtained by fitting all data points (i.e., HI bins) from Fig. lAll The errors on the parameters of the fit were determined by bootstrap resampling the 
quasars. Dotted curves indicate the Icr confidence limits. 



ear least-squares fit to the data gives a = 0.71^Q [7+0.56^o 4g x 
(z-3) (where the errors were again computed by bootstrap re- 
sampling the quasars). As expected, the two methods of mea- 
suring the gradient of the metallicity with density give very 



similar results. 

Let us now turn to the width of the lognormal fit to the 
metallicity distribution. Figure [121 shows er([C/H]) versus 
S (left-hand panel) and versus z (right-hand panel). The 
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FlG. 1 1 . — Median metallicity vs. overdensity 8 for different redshift bins. From the top l eft t o the bottom right panels correspond to the following redshift 
ranges: z = 1.5-2.5, 2.5-3.5, and 3.5-4.5. Data points with lcr errors were taken from Fig. lAll Light-colored data points have a lower lcr error of -co. Solid 
lines indicate least-squares fits to the data points. Dotted curves indicate the lcr confidence limits, which were computed by bootstrap resampling the quasars (for 
the z = 4 bin, which contains only two quasars, we resampled the data points). 



data points correspon d dire ctly to the points conn ected by the 
dashed lines in Figure lATl As discussed in 35.21 our method 
for computing cr([C/H]) relies on fitting different percentiles 
of tciv for a fixed thi, and because the percentiles are cor- 
related, the derived error bars on er([C/H]) may not be ro- 
bust. However, the relative errors among the different points 
should still be reliable and hence the least-squares fits to the 
data points — shown as the solid lines in Figure^] — should 
be accurate. The errors on the fits were computed by bootstrap 
resampling the quasars and should therefore also be robust es- 
timates of the statistical errors. Note, however, that the fact 
that the measurements of the scatter are mostly one-sided (the 
low-metallicity tail is undetected), could result in somewhat 
larger systematic errors. 

We find a significant (> 2a) decrease in the scatter 
with overdensity: ct([C/H]) = 0.76^;g|-0.24^;^(log (5-0.5). 
There is a hint of evolution in the data: the best-fit scatter de- 
creases by about 0. 14 dex from z = 4 to 2, but the a constant 
value is acceptable at about the 1.2er level. 

The anticorrelation of er([C/H]) with log S is somewhat dif- 
ficult to see by eye, because the trend is not large compared 
with the scatter between the data points. There are, however, 
enough points to detect weak anticorrelations. The Spearman 
rank-order correlation coefficient is -0.24, indicating that the 



anticorrelation exists at 97% confidence. Note that this test is 
based on a ranking of the data points and does not make use 
of the errors. 

We can summarize all of our data by fitting functions of red- 
shift and overdensity to the median metallicity and the (log- 
normal) scatt er, using all the log<5 > -0.5 data points shown 
in Figure lATI The function 

[C/H]=-3.47!^ + 0.08t^ x (z-3) 

+0.65!^ x (log (5 -0.5) (8) 

provides a very good fit to the data: \ 2 = 1 14. 1 for 1 84 degrees 
of freedom. As for the projection fits, the reduced \ 2 is l° w 
because we overestimated the errors for the data points with 
tciv ~ t c , i.e., points for which C IV is barely detected (see 
^5.11 . Excluding all data points with thi < t c does not change 
the fit (the differences are smaller than the statistical errors), 
but does result in a much more reasonable reduced \ 2 of 0.89. 

To test whether there is evidence for density-dependent 
evolution or, equivalently, a redshift-dependent gradient with 
density, we have also fitted a function including a (z — 
3)(log<5-0.5) term to the data. The best-fit coefficient of the 
cross term is consistent with zero at the w 0.8<7 level, and 
including it improves the quality of the fit only slightly, to 
X 2 = 112.5 for 183 dof. 
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FIG. 12. — Standard deviation of th e log normal metallicity distribution as a function of overdensity (left) and redshift (right). Data points with 1<t errors 
correspond to the HI bins shown in Fig. lAll The relative sizes of the errors are reliable, but the absolute sizes are not. Solid lines are least-squares fits to the data 
points. Dotted curves indicate the la confidence limits, which were computed by bootstrap resampling the quasars. 



Table 2. Fits for [C/H] = a+/3(z-3)+7(log<5-0.5) 



» r 11 +ler.2(7 r,+l<T,2er +\a,2a , r 

Model Q_ 1ct , 2ct P-la.la l-ia,2a A.O.I. 

C\r 1 A7+007.0.13 nns tO.fl9.0.1! „ ,-+0.10,0.18 

QG - 3 - 47 -0.06.0.12 + 008 -0.1 0.0.20 +0 - 65 -0.I4,0.27 114.1,184 

„ , q ,+0.07,0.I3 nn ,+0.09.0.21 „ ,.,+0.08,0.15 .,,„ , s4 
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nr5 n , 7e +0.14,0.21 n , -+0.06,0.13 , n Q ,,+0.I8,0.28 , 7 

QGS3.2 -3.78_ 09 „ 16 -OM^ n 022 +0.93_ 32 47 1 24.2, 184 

nrs .,.+0.06,0.11 „ +0.1 0.0.21 .-.+0.07,0.17 , „ . 

QGS -4- 1 4 - .05. a u +Q - 54 -o.07.o.i6 +L31 -o.07.o.i8 114.2,184 



Table 3. Fits for a = Q+-/3(z-3)+7(log<5-0.5) 



Model Q , + 1(T ' 2 < T o+Ict.2,7 -Her, 

r>r- , mt* 05 ' ' 09 , n m+0.08.0.17 n OQ +0.( 

QG -+0-76-0.08.0.16 +002 -o.i2.o.22 -°- 23 -o.( 
Q As for QG 

QGS3.2 As for QG 

QGS As for QG 



Fitting a three-parameter surface to the scatter in the log- 
normal distribution yields 

<r([C/H]) = 0.763$ + 0.02!ft°t x (z-3) 

-0.233$ x (log 6- 0.5). (9) 

As was the case for the median, using the four-parameter 
function does not improve the fit and yields no evidence for a 
nonzero coefficient of the (z-3)(log 6-0.5) term. 

The best-fit parameters for all the surface fits are listed to- 
gether with their 1 and 2er errors in Tables|5]and|5] Note that 
all these fits are based on data in the range log<5 = -0.5- 1.8 
(log* = -0.5-0.9 for z> 3.5), z= 1.8-4.1 (z = 2.1 -3.3 for 
log 5 > 1) and that extrapolations outside this range of param- 
eter space could be inaccurate. 

The surface fits confirm the picture suggested by the projec- 
tion plots. There is very little room for evolution, but there is 
strong evidence for both an increase in the median metallicity 
and a decrease in the lognormal scatter with overdensity. 

7.1. Varying the UV background 



Al l ionization corrections were computed as discussed in 
jQI using model QG (Haardt & Madau 2001) for the UV/X- 
ray background from galaxies and quasars (see 34. 2> . rescaled 
to match the evolution of the mean Lya absorption as mea- 
sured from our sample of observations (see {j3}. Since the 
spectral shape of the UV background is not well constrained, 
it is important to investigate what effect changes in the spec- 
trum have on the derived metallicities. 

In particular, since photons with energies greater than 4 ryd 
can ionize C IV but not C III, the ionization corrections can 
be sensitive to the break at the He II Lyman limit, which 
is much greater if the background is dominated by stars as 
opposed to quasars, and can also be very large if He II is 
not fully reionized. As Figure [3] shows, for gas with densi- 
ties «h ^ 10~ 5 cm" 3 the z = 3 ionization correction is much 
larger for model Q — which includes only contributions from 
quasars — than for model QG. Thus, as expected, for low gas 
densities the derived metallicity is sensitive to the hardness 
of the UV background radiation, with harder spectra yielding 
higher metallicities. 

To investigate how sensitive the metallicities are to changes 
in the UV radiation field, we have recomputed the surface fit 
[C/H](S,z) (eq. 0) using the two add ition al models for the 
UV background radiation described in 34.21 Model Q (Haardt 
& Madau 2001) is an updated version of the Haardt & Madau 
(1996) model for the background radiation from quasars only. 
Model QGS is identical to model QG, except that above 4 ryd 
the flux has been decreased by a factor of 10. The results are 
listed in Table [2] 

Model QGS is much softer than model QG and therefore 
gives lower metallicities and a much stronger gradient with 
overdensity. In addition, the metallicity is everywhere de- 
creasing with time, clearly an unphysical 18 result. We get 
the same unphysical trend if we decrease the flux above 4 ryd 
by only a factor of 3 instead of a factor of 10 as in QGS. 
Although using model QGS at all redshifts results in a metal- 
licity evolution that does not make physical sense, QGS could 
be a reasonable model for the UV background at z ~ 4 if He II 

18 Note that for fixed high overdensities the metallicity c ould in principle 
decrease with time because of infall of metal-poor gas (see t|8.1l . 
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was reionized late. We have therefore tested model QGS3.2, 
which is identical to model QG for z < 3.2 and identical to 
model QGS for z > 3.2. This transition redshift was chosen 
because various lines of evidence suggest that He II reion- 
ization was incomplete at earlier times (e.g., Songaila 1998; 
Schaye et al. 2000b; Heap et al. 2000; Theuns et al. 2002a). 

Figure^]shows the two projections of the three-parameter 
surface fits (listed in Table|2) for models QG, Q, and QGS3.2. 
The left-hand panel shows the median metallicity versus log S 
for z = 3, while the right-hand panel plots the median metal- 
licity versus z for log<5 = 0.5. Although these plots cannot 
shed light on density-dependent evolution, they do illustrate 
the key differences between the various UV backgrounds. 

As for model QG, model Q gives no evidence for evolution. 
However, for model Q the results do differ from our fiducial 
model in other ways. Model Q yields much higher metallic- 
ities for log<5 < 1 and gives a much smaller gradient of the 
median metallicity with density. In the next section we will 
show that this small gradient, combined with the observed de- 
crease in the scatter with density, results in a mean metallicity 
that decreases significantly with increasing density, probably 
an unphysical result. Note that model Q also does not fare as 
well as model QG in fitting the observed C III/C IV ratios (see 

Model QGS3.2 yields a significantly lower metallicity for 
z > 3.2 than our fiducial QG model. As a result, there is ev- 
idence for (weak) evolution (d[C/U]/dz = -0.12^;?f). Be- 
cause many of our low-density detections have high redshifts, 
the gradient with overdensity is even stronger than for model 
QG. Since it is only the high-redshift data that are responsible 
for the larger gradient, there is evidence for density-dependent 
evolution, including a cross term in the fit improves it by 
A% 2 = 9 relative to the three-parameter fit. 

Although uncertainties in the UV background are important 
for estimates of the median metallicity, this is not the case for 
our measurements of the lognormal scatter, whic h are inde- 
pendent of the assumed UV background (see i5.2\ . Our find- 
ing that the metallicity distribution is lognormal and our mea- 
surements of the scatter in this distribution are therefore much 
less model-dependent than our measurements of the median 
metallicity. They are, however, still not completely model- 
independent because we have assumed that the UV radiation 
is constant for fixed thi. If there were significant scatter in the 
UV radiation field corresponding to a fixed overdensity, then 
this would generate scatter in the rciv distribution, leading us 
to overestimate the width of the metallicity distribution. 

8. DISCUSSION 
8.1. Summary of measurements 

We have measured the distribution of carbon as a function 
of overdensity and redshift using data in the range log 5 = 
-0.5- 1.8 and z= 1.8-4.1. For a fixed overdensity and red- 
shift the metallicity distribution is close to lognormal, at least 
from about -0.5 to +2a. We measure a lognormal scatter 
a([C/H]) = 0.76!° : °i + 0.02^?|(z-3)-0.23^(log5-0.5). 
Thus, we find no evidence for evolution, but we do find a 
significant decrease in the scatter with overdensity. The mea- 
surements of the scatter are independent of the spectral shape 
of the UV background radiation. 

Unlike the scatter, the median metallicity does depend on 
the model for the UV background, although it is insensitive 
to the spectral shape for S > 10. For our fiducial model QG, 
which includes contributions from both galaxies and quasars 



(see Haardt & Madau 2001), we find [C/H] = -3.47^ + 
0.08^(z-3) + 0.65!j] ^(log(5-0.5), i.e., no evidence for 
evolution, but a strong gradient with overdensity. 

Harder UV backgrounds yield higher metallicities at low 
overdensities. Model Q, which includes only UV radiation 
from quasars, gives a median metallicity higher by about 
0.6 dex for log<5 = 0.5, a much weaker gradient with over- 
density (d[C/E]/dlogS = 0.17^|), but again no evidence 
for evolution. As we will discuss below, this UV back- 
ground is probably too hard since it yields a mean metallicity 
that decreases with overdensity. Note that measurements of 
Nudi/Nm also imply that the true UV background is softer 
than model Q (Heap et al. 2000; Kriss et al. 2001). Mak- 
ing the UV background too soft also gives unphysical results: 
model QGS, which has a 10 times smaller flux above 4 ryd 
than model QG, 19 yields a metallicity that strongly decreases 
with time for all overdensities. We find the same unphysi- 
cal trend if we decrease the flux above 4 ryd by a factor of 3 
instead of 10. 

Hence, it appears that the metal distribution (median and 
scatter) evolves very little from z = 4 to z = 2, which sug- 
gests that most of the enrichment of the low-density IGM took 
place at higher redshifts. There is, however, a caveat: we do 
find (weak) evolution in the median metallicity if we make the 
UV background much softer at high redshift (or much harder 
at low redshift). For example, model QGS3.2, which has a 
10 times smaller flux above 4 ryd for z > 3.2 (as may be ap- 
propriate if He II had not yet reionized by that time), gives 
d[C/R]/dz = -0.12^ for log<5 = 0.5 and stronger evolution 
for lower overdensities. 

Note that even if the enrichment of the IGM was completed 
at z ^> 4, one would still expect some evolution in the gradient 
of metallicity with overdensity, unless the initial gradient was 
zero. Because gravity tends to increase density contrasts, any 
initial gradient of metallicity with overdensity should become 
weaker with time. In fact, our measurements favor a small 
increase in the gradient with redshift: for model QG the co- 
efficient of the (z— 3)(log 5-0.5) term in our four-parameter 
surface fit is 0.23^. Although a detailed comparison with 
models for the enrichment of the IGM is beyond the scope of 
this paper, we note that the simulations we have performed 
show that the expected "passive" evolution in the gradient of 
the metallicity with overdensity from z = 4 to z = 2 is suffi- 
ciently small to be easily compatible with our measurements. 

8.2. Mean metallicities 

Having measured the median and the scatter of the lognor- 
mal distribution, we can also compute the mean (see eq. Q). 
For our fiducial model we find log (Z) -2.8 at log 5 = 0.5, 
where Z is short for 10 [C / HI , with no evidence for evolution. 
At the 1 a (2a) level the data allow an increase of about 0. 1 dex 
(0.3 dex) per unit redshift. For model Q, log (Z) is ~ 0.6 dex 
higher at this overdensity and, as for model QG, there is no 
evidence for evolution. Note that the limits on the allowed 
evolution of the mean metallicity are weaker than for the me- 
dian metallicity because there is additional uncertainty in the 
redshift dependence of the lognormal scatter. 

For our fiducial model the mean metallicity increases with 
overdensity from log (Z) = -2.9 at 5 ~ 1 to -2.3 at 5 ~ 10 2 for 
z = 3. For model Q, however, it decreases from -2.1 to -2.5 
over the same density range. This negative gradient is signifi- 

19 Reducing the flux above 4 ryd further has very little effect because the 
CIV ionization rate is already small. 
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FIG . 13. — Visualization of the fits to the median metallicity as a function of overdensity and redshift for z = 3 (left-hand panel) and log <5 = 0.5 (right-hand panel) 
for our fiducial UV background model QG, as well as for models Q and QGS3.2. The shaded regions enclose the Icr confidence contours for the three-parameter 
surface fits listed in Tablel2l Harder UV backgrounds yield higher metallicities for S < 10. Strong evolution is ruled out for all models. 



cant at the 1 a (but not 2a) level and may imply that model Q is 
unphysical since models for the enrichment of the IGM gener- 
ically predict positive metallicity gradients with overdensity 
(e.g., Cen & Ostriker 1999; Aguirre et al. 2001a, 2001b). 

8.3. Global metallicities 

By combining our measurements of the carbon distribu- 
tion (equations JS) and |9)) with the mass-weighted probabil- 
ity density distribution for the gas density obtained from our 
hydrodynamical simulation, we can compute the contribution 
of the carbon that we see to the global, mass-weighted mean 
metallicity. For z = 3 and log<5 = -0.5-2.0 we find a cosmic 
carbon abundance [C/H] =-2.80 ±0.1 3, with no evidence for 
evolution. Extrapolating our measurements to the full density 
range of the simulation, we find [C/H] = -2.56 ±0.16 (the 
higher overdensities being responsible for the difference). For 
model Q all values are about 0.5 dex higher. 

Relative to the critical density, our measurement of the con- 
tribution of the forest to the global carbon abundance corre- 
sponds to a global carbon density of 

Oc=10 [C/H] + (C/H) Omc «H ; (10) 
Pc 

«2.3xl0- 7 10»' 8 (^), (11) 

where ;«c is the atomic mass of carbon, «h is the total, comov- 
ing number density of hydrogen, and p c is the critical mass 
density at redshift zero. 

8.4. Filling factors 

The dashed curves in Figure [141 show the fraction of gas 
with a metallicity greater than [C/H] as a function of [C/H] 
for various overdensities. All curves are for model QG and 
z = 3. The filling factors were computed using the surface fits 
to the median metallicity and lognormal scatter as a function 
of overdensity and redshift (equations |8) and 0). Note that 
for metallicities below the median (i.e., filling factor 0.5), we 
detect C IV only for significantly overdense gas. Thus, the 
lognormal shape of the curves in the upper half of the plot 
is a result of extrapolating the shape measured for the bot- 
tom half. The thick, solid curve shows the total filling factor, 
i.e., the volume fraction. It was computed by combining our 



metallicity measurements with the volume weighted proba- 
bility density distribution for the gas density in our hydrody- 
namical simulation at z = 3 (using a physical smoothing scale 
of 75 kpc). The filling factors for z = 2 and z = 4 (not plotted) 
fall nearl y on top of this curve. 

Figure [141 summarizes many of our results. The intersec- 
tion of the dashed curves with the filling factor 0.5 gives the 
median metallicity at the various overdensities. The fact that 
the dashed curves become steeper going from low to high 
overdensities {left to right) reflects our finding that the scatter 
decreases with overdensity. One consequence of this is that 
high-metallicity gas ([C/H] ^ -2) is rare 20 . For all metallici- 
ties the filling factor increases with overdensity. For model Q, 
which gives a median metallicity that is nearly independent of 
overdensity, this is in fact not the case for high metallicities. 
Since a filling factor that decreases with overdensity seems 
implausible, this again suggests that the UV radiation is too 
hard in model Q. 

Essentially all collapsed gas clouds (5 > 10 2 ) have a metal- 
licity greater than 5 x 10 -4 Z Q . Metallicities of 1 and 5 times 
10~ 4 have been claimed to be the maximum possible metal- 
licities for the formation of supermassive stars (Bromm et al. 
2001; Schneider et al. 2002). Our results therefore indicate 
that by z = 4 the formation of such stars had already come to 
an end. Some caution is appropriate, however, because the top 
parts of the dashed curves are based on extrapolations since 
we generally detect C IV only if the metallicity is similar or 
greater than the median value. 

8.5. Do we see all of the carbon? 

Throughout this paper we assumed that the C IV gas has a 
temperature T <C 10 5 K, so that collisional ionization is unim- 
portant. Our measurements of the C III/C IV ratio strongly 
support this assumption for log S > 1 , but cannot rule out high 
temperatures for lower density gas. However, the fact that we 
detect C IV generally even for H I optical depths as low as 
thi ~ 1 implies that hot gas would need to have a substantial 
filling factor to invalidate the assumption that the gas that we 
observe is predominantly photoionized. Note that Carswell et 
al. (2002) and Bergeron et al. (2002) could rule out collisional 

20 There could be somewhat more high metallicity gas than suggested by 
this figure becaus e we do not have enough pixels to probe the metallicity 
beyond +2a (c.f. i|5.2l . 
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FIG. 14. — Filling factor of gas with a carbon abundance greater than [C/H] 
plotted as a function of [C/H] for z = 3 and UV background QG. Dashed 
curves indicate the filling factors at overdensities of, from left to right, log S = 
-0.5,0.0, . . . ,2.0. The thick solid curve indicates the total filling factor, i.e., 
the volume fraction. 



ionization for many of the O VI absorbers associated with low 
column density H I systems on the basis of their line widths. 

Gas with a temperature T 3> 10 5 K is too highly ionized to 
cause detectable absorption in C IV or H I and will therefore 
not introduce errors in our measurements of the metallicity of 
the warm, photoionized component. It is, however, important 
to note that our analysis cannot reveal carbon hidden in hot, 
X-ray-emitting gas, which may contain a large fraction of the 
intergalactic metals (Theuns et al. 2002b). The same is true 
for cold, self-shielded gas for which most of the carbon is 
neutral or singly ionized. Comparisons to models in which a 
significant fraction of intergalactic carbon is in gas with T ^> 
1C I 5 K or T <C 10 4 K must therefore be made with care. 

8.6. Comparison to previous work 

There have been no previous systematic attempts to mea- 
sure the metal distribution in the Lya forest as a function of 
overdensity and redshift. Previous studies using C IV pixel 
statistics (Cowie & Songaila 1998; Ellison et al. 1999, 2000) 
reported optical depth measurements consistent with our re- 
sults. However, these studies used less sensitive methods, 
considered only the median Tcrv> an d did not attempt to cor- 
rect for ionization. Published measurements of the carbon 
abundance have generally been based on a comparison of 
Voigt profile statistics of C IV lines with either a hydrody- 
namical simulation assuming a uniform metallicity or simple 
photoionization models for individual absorbers. 

Both simulations (e.g., Haehnelt et al. 1996, Rauch et al. 
1997; Hellsten et al. 1997; Dave et al. 1998) and simple 
photoionization models (e.g., Cowie et al. 1995; Songaila 
& Cowie 1996; Rollinde et al. 2001; Carswell et al. 2002) 
find that a carbon abundance of 10~ 3 to 10~ 2 solar repro- 
duces the z ~ 3 observations. Rauch et al. (1997), Hellsten 
et al. (1997), and Dave et al. (1998) all find [C/H] = -2.5 
and also report evidence for scatter of about 1.0 dex (Dave 
et al. 1998 find 0.5 dex). These studies used the relatively 
hard Haardt & Madau (1996) model for the UV background 
from quasars, which closely resembles our model Q. For 
this model we find a median metallicity at z = 3 of [C /H] = 
-2.91^ + 0.17^;g|(log 5-0.5) and a lognormal scatter of 



cr([C/H]) = 0.76^ -0.23^$(log 6 -0.5), both nearly inde- 
pendent of redshift. However, compared with our work, the 
metallicity measurements of previous studies were dominated 
by denser gas: log (5 ~ 1.0- 1.5. For these densities we find a 
scatter of about 0.6 dex, a median metallicity [C/H] w -2.8, 
and a mean metallicity log ([C/H]) « -2.4 (see eq. Q; note 
that it is unclear whether we should be comparing the mean 
or the median when comparing with previous work). Hence, 
if we use a hard UV background and focus on the high den- 
sity gas (Nm ~ 10 15 - 10 16 cm" 2 ), then our results for z = 3 
agree with the metallicity measurements obtained by previous 
studies. However, our fiducial UV background (QG) gives a 
smaller metallicity, particularly at low overdensities. 

Songaila (2001) measured the C IV column density distri- 
bution as a function of redshift from z = 1.5 toz = 5.5 and 
found no evidence for evolution. It is, however, important not 
to overinterpret this finding. To draw conclusions about the 
evolution of the metallicity, one must correct for ionization, 
and this correction is time-dependent because the universe ex- 
pands and the UV background evolves. 

Songaila (2001; see also Pettini et al. 2003) also computed 
ficiv as a function of redshift by summing all of the observed 
C IV systems. This provides a strict lower limit to be- 
cause not all carbon is triply ionized and because some C IV 
systems may have been lost in the noise and contamination. 
Finite sample size is also important because the integrals of 
fits to the C IV column density distribution diverge at high 
N C iv- Songaila found fi C iv(10 12 cm -2 < New < 10 15 cm -2 ) w 
(3 ±2) x 10~ 8 at all redshifts. 21 Using our measurements of 
the carbon distribution for model QG and the mass-weighted 
probability density distribution for the gas density from our 
hydrodynamical simulation, we find £Iq ~ (2.3 ±0.7) x 10~ 7 
for \o gS = -0.5-2.0, with no evidence for evolution (see 
38. 3> . For model Q the cosmic carbon abundance is about a 
factor of 3 higher, again with no evidence for evolution. Thus, 
we find that fic is indeed much higher than ficiv measured 
by Songaila (2001). The large difference implies that the ob- 
served evolution of ilciv tells us very little about the evolution 
of^ c . 

Nevertheless, Songaila's measurements of Slav provide an 
important consistency check. Using our measurements of the 
distribution of carbon, the density distribution of our hydro- 
dynamical simulation, and our ionization correction as a func- 
tion of overdensity and redshift, we can compute flciv- For 
-0.5 < log S < 2.0 we find Ociv/lO" 8 « 1.3 ±0.8, 4± 1, and 
8 ± 5 for z = 2, 3, and 4, respectively (as expected, the differ- 
ent UV background models give similar numbers). The values 
are nearly the same if we extrapolate to higher or lower over- 
densities, which implies that we are seeing most of the C IV. 
The quoted errors are estimates of the statistical errors; i.e., 
they are based on the errors on the parameters of the surface 
fits. Note that the difference between ficiv and J7c decreases 
from z = 2 to z = 4. Thus, our results are in good agreement 
with those of Songaila (2001). Given that our calculation of 
flew relies on convolving our measured carbon distribution 
with the gas distribution extracted from our simulation, this is 
a highly nontrivial consistency check. 

21 Songaila (2001) measured C C 1V ~ (5 ±4) X 10~ 8 for z = 2-4 assum- 
ing Q m = 1 . Since her method is based on measuring the total CIV column 
density per unit redshift, this value scales as \J Q, m + (1 — f2,„)/(l +z) 3 for a 
cosmologically flat universe. For our cosmology (Q„, = 0.3) this becomes 
ficiv ~ (3 ±2) X 10~ 8 . 



Distribution of Metals as Traced by C IV 



19 



8.7. Uncertainties 

Our method contains several steps and parameters chosen to 
minimize the errors in the metallicity measurements in simu- 
lated spectra. To test the robustness of our method to these 
choices, we have redone the surface fits after dividing each 
quasar spectrum in two, doubling the H I bin size, omitting 
our correction to the continuum fit of the C IV region, exclud- 
ing all data with 5 < 1, excluding all data with thi < r c , or 
excluding all H I bins containing fewer than 50 pixels and/or 
contributions from fewer than 10 different chunks (our fidu- 
cial values for these parameters are 25 and 5, doubling them 
removes many of the higher H I bins). We find that the co- 
efficients of the various surface fits always agree within their 
errors with those obtained using our standard method. Thus, 
our results are insensitive to the details of our method and the 
surface fits are not determined by the highest/lowest H I bins. 

Although our results are robust with respect to small 
changes in the methodology, there could of course still be 
systematic errors. The main uncertainty in the median metal- 
licity comes undoubtedly from the uncertainties in the spec- 
tral shape of the UV background. Although the lognormal 
scatter is independent of the mean spectral shape, it is sen- 
sitive to fluctuations. If there were significant fluctuations in 
the UV background, then this would introduce scatter in the 
tciv(thi) relation even for a uniform metallicity, particularly 
at low overdensities where the ionization correction is large. 
If fluctuations are important, then our measurements should 
be interpreted as upper limits on the scatter in the metallicity. 

Even for a uniform UV background there are likely to be 
some systematic errors in our ionization corrections. For ex- 
ample, our measurements of the C III/C IV ratio cannot rule 
out the possibility that collisional ionization is important for 
6 <C 10. At some level, the ionization corrections must be in- 
correct, since redshift-space distortions will result in multiple 
gas elements contributing to the absorption at a fixed redshift. 

Fortunately, we can use our hydrodynamical simulation to 
estimate the size of these and other systematic uncertainties: 
we can distribute metals in our simulation according to our 
measurements, create synthetic versions of our observations, 
and then co mpare the simulated and observed optical depth 
statistics. In 35. 21 we showed that our method passed this test 
for Q1422+230. It is important to note that interpreting the 
results of this test is not completely straightforward, because 
there is some ambiguity in the interpretation of our measure- 
ments. For example, we do not know the exact scale on which 
our o verdensities are smoothed (probably ~ 50- 100 kpc; see 
34. 3L and we do not know the exact scale(s) over which we 
are measuring the scatter. Consequently, it is not completely 
clear how to distribute the metals in the simulations. Having 
said this, we expect the density smoothing scale to be similar 
to the smoothing scale of the forest — which depends on the 
gas density and is not very different from the resolution of our 
simulation — and we find that the precise scale on which the 
scatter is imposed is not critical provided it is greater than the 
smoothing scale of the forest. 

We ha ve ca rried out the test described above (and in more 
detail in 35.21 for our entire sample, distributing the metals in 
our simulations according to equations (|8j and l|9}, although 
we neglected the small, observed decrease in the scatter with 
overdensity since we have no practical way of implementing 
such a trend in the simulation (we imposed the scatter on a 
scale of 1.2/i _1 Mpc comoving, i.e., l/10th of the simulation 
box size). Averaging over 10 realizations, we find for the data 



points with thi > r c a reduced \ 2 of 1.02, 0.97, and 1.35 for 
the 50th, 69th, and 84th percentiles, respectively. The total 
number of data points included in this estimate is 92 for the 
50th percentile and 95 for the other two percentiles. Hence, 
we have a model that can reproduce the observed optical depth 
statistics. It appears that any systematic errors must be small 
and that the ambiguities discussed above are not critical. 

We also carried out the following, more stringent test. As 
before, we first generate a set of synthetic observations us- 
ing our three-parameter fits to the data (eqs. |8) and (9), but 
with da/dlogS = 0). Second, we measure the carbon abun- 
dance from the simulated spectra in the same way as we did 
for the data, and we fit the same three-parameter functions to 
the results. Finally, we compare the fits to the ones we started 
with. The differences between the fits give us estimates of 
the systematic errors, or rather upper limits to the systematic 
errors given the ambiguities in the smoothing scale discussed 
above. Averaging over 10 realizations we find systematic er- 
rors for the coefficients of the median metallicity (eq. 0) 
of Aa = -0.10 (1.7ct), A/3 = -0.20 (2.0cr), and A<7 = -0.14 
(I.Oct), where the numbers in brackets give the ratios of the 
systematic to the statistical errors. The variance between the 
different realizations is comparable to the statistical errors. 
Note that the signs of the systematic errors is such that we may 
have slightly overestimated the evolution (i.e., d[C/H]/dz too 
negative) and underestimated the gradient with overdensity. 
For the coefficients of the scatter (eq. |9)) we find Aa = -0. 1 1 
(0.69er), A/3 = 0.02 (0.25<r), and A 7 = 0.01 (0.1 la). Thus, 
the systematic errors are small for the scatter, but may not 
be insignificant for the normalization (a) and the trend with 
redshift (/3) of the median. No systematic errors are greater 
than twice the corresponding statisical errors. Although these 
tests give us confidence that our results are robust, they do of 
course not rule out the existence of other, unknown, sources 
of systematic errors. 

8.8. Future prospects 

How could our measurements be improved? Higher signal- 
to-noise ratios (and better continuum fits, but these are signal- 
to-noise limited) are required to measure the median metallic- 
ity at lower overdensities and to verify that the metal distribu- 
tion remains lognormal below the median metallicity. Better 
coverage of the higher-order Lyman lines would enable us to 
increase the number of data points for log 8 > 1 . To improve 
the constraints on the evolution more data at z > 3 would be 
very helpful since only four out of our 19 quasar spectra have 
a median redshift greater than 3. For S < 10 the median metal- 
licity becomes increasingly sensitive to the softness of the UV 
radiation, and it will be necessary to improve the constraints 
on the spectral shape of the UV background to tighten the 
constraints. 

9. CONCLUSIONS 

We have measured the statistical correlation between C IV 
and H I absorption in a sample of 19 high resolution 
VLT/UVES and Keck/HIRES spectra of QSOs spanning a 
redshift range ~ 1.5-4.5. Hydrodynamical simulations 
(which can reproduce the observed statistics of the Lya for- 
est) link the observed H I absorption to the density and tem- 
perature of the absorbing gas. Given a model of the ionizing 
background radiation field, this allows us to convert the me- 
dian ratio of C IV to H I pixel optical depths into a median car- 
bon abundance [C /H] as a function of density and redshift. In 
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Paper I we tested this technique extensively on realistic sim- 
ulated spectra, and here we have generalized the method to 
recover the distribution of metals from the full distribution of 
tciv as a function of r H i. 

Our fiducial model for the UV background is that of Haardt 
& Madau (2001) — which includes contributions from both 
galaxies and quasars — renormalized using our measure- 
ments of the evolution of the mean Lya absorption. 

Our primary conclusions are as follows: 

• The best-fit power-law evolution of the effective Lya 
optical depth in the spectra of 21 quasars is logr e ff = 
(-0.44 ± 0.01) + (3.57 ± 0.20)log[(l +z)/4]. Contam- 
ination from metal lines — which we removed before 
making these measurements — is significant, particu- 
larly at lower redshifts. 

• C III to C IV optical depth ratios typical of collision- 
ally ionized C IV (T > 10 5 K) are ruled out by the 
data, at least for overdensities 6 > 10. Simulations 
with T r~j 10 4 K naturally reproduce the observations. 
The existence of very hot (T 3> 10 5 K) or very cold 
(r <C 10 4 K) carbon cannot be ruled out as it would not 
give rise to significant C III or C IV absorption. 

• For a fixed overdensity 6 and redshift z the metallicity 
distribution is close to lognormal, at least from about 
-0.5 to +2a. 

• Over the range log S = -0.5 - 1 .8 and z = 1.8-4.1 the 
median carbon abundance is well fitted by the fol- 
lowing function: [C/H]=-3.47!g;^ + 0.08!g$(z-3) + 
0.65^0 ["(log 5-0.5). The overdensities have effectively 
been smoothed on the same scale as the Lya forest 
(<~ 50- 100 kpc). From tests using simulations we find 
that the known systematic errors are at most twice as 
large as the quoted statistical errors. There is no evi- 
dence for a nonzero (z- 3)(log<5-0.5) term. Thus, we 
find no evidence for evolution, but strong evidence for 
a positive gradient with overdensity. 

• Harder UV backgrounds give higher carbon abun- 
dances for log S < 1 and thus smaller gradients with 
overdensity. The (rescaled) Haardt & Madau (2001) 
quasar only UV background yields a mean metallic- 
ity that decreases with overdensity, which is probably 
unphysical. UV backgrounds that are much softer than 
our fiducial model result in a median (and mean) metal- 
licity that strongly increases with redshift, which is also 
unphysical. However, a UV background that is softer 
only for z > 3.2, as may be appropriate if the reioniza- 
tion of He II was incomplete before this time, is allowed 
and results in significant (but still weak) evolution in the 
median metallicity. 

• We find a nonzero median carbon abundance for un- 
derdense gas (log 6 = -0.5-0.0) at the 99.2% (2.4a) 
confidence level. 

• The lognormal scatter is ct([C/H]) = 0.76!g^ + 
0.02!™f (z - 3) - 0.23!g$(log£ - 0.5). Our measure- 
ments of the scatter are strictly independent of the 



spectral shape of the UV background radiation. They 
are, however, sensitive to fluctuations in the UV back- 
ground, which we have assumed to be small. Tests us- 
ing simulations show that the (known) systematic errors 
are smaller than the quoted statistical errors. Thus, we 
find no evidence for evolution, but we do find strong 
evidence for a modest decrease in the scatter with over- 
density. 



• By combining our measurements of the distribution of 
metals with the mass weighted probability density dis- 
tribution in our hydrodynamical simulation, we find 
that the gas between log<5 = -0.5 and 2.0 accounts 
for a cosmic carbon abundance of [C/H] = -2.80 ± 
0.13, with no evidence for evolution. Expressed in 
terms of the critical density, this becomes ~ 2 x 

10" 7 (£V0.045). We find that n c > ^civ which im- 
plies that ionization corrections are important. There- 
fore, the finding that Slav is roughly constant between 
z = 2 and z = 5 (Songaila 2001) does not by itself pro- 
vide interesting constraints on the evolution of fie- 



We have presented the first measurements of the distri- 
bution of carbon as a function of overdensity and redshift. 
Published models for the enrichment of the IGM by galac- 
tic winds and/or Population III stars predict a wide variety of 
metal distributions (e.g., Cen & Ostriker 1999; Aguirre et al. 
2001a, 2001b; Madau, Ferrara, & Rees 2001; Theuns et al. 
2002b; Scannapieco, Ferrara, & Madau 2002; Thacker, Scan- 
napieco, & Davis 2002; Furlanetto & Loeb 2002; Scanna- 
pieco, Schneider, & Ferrara 2003), which is perhaps not sur- 
prising given the large uncertainties in the parameters of the 
models. Our measurements will provide some much needed 
constraints and it will be interesting to see how they narrow 
the parameter space for the models. 

More generally, the measurements reported here should 
provide a new and stringent set of targets for the next gen- 
eration of models for the history of the formation of stars 
and galaxies, and the effects of feedback on the intergalactic 
medium. 
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APPENDIX 
A. RESULTS FROM INDIVIDUAL QUASARS 

Figure IaTI shows the results for all 19 quasars listed in Tabled For each quasar two panels are plotted. The left-hand 
panels show the median recovered C IV optical depth as a function of the recovered thi, while the right-hand panels show the 
median metallicity as a function of overdensity. 22 The dashed lines in the right-hand panels indicate the ±1ct scatter measured 
from lognormal fits to the distribution of C IV pixel optical depths. Our methods for measuring the median and variance of the 
metallicity distribution are explained in |5] 

To facilitate future compari sons with models for the distribution of carbon, we have tabulated the optical depths for several 
computed percentiles in Table lAll 



22 Recall that all overdensities have effectively been smoothed on the same scale as the Lya forest spectra. As discussed in SI4.3I the smoothing is due to a 
combination of thermal broadening and Jeans smoothing and varies from ~ 30 to 100 kpc (physical) depending on density and redshift. 
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FIG. A 1 . — Results for individual quasars. For each quasar two panels are plotted. Left: Median recovered CIV optical depth as a function of the recovered 
HI optical depth. The horizontal dotted line indicates r mm ; the \u error in this level is indicated on the right-hand side of the panel. The solid line is th e best-fit 
power law, eq. |2}- Right: Median carbon abundance [C/H] as a function of log overdensity. All metallicities were computed as described in j|5.1l using the 
QG UV background. The best-fit power law to the Tffl > t c data points is shown as the solid line (dotted where extrapolated). The dashed lines show the ±1<t 
lognormal scatter as measured from the full distribution of pixel optical depths (see i!5.2l . All data points are shown with both 1 and la error bars. 
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Table B 1 . Observed effective optical depth 
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±0 


037 


Q1425+604 


2 


694 


2 


544 


2 


844 





417 


±0 


044 


2 


670 


2 


544 


2 


882 





232 


±0 


023 


Q1425+604 


2 


994 


2 


844 


3 


144 





363 


±0 


034 


3 


009 


2 


882 


3 


144 





360 


±0 


040 


Q2126-158 


2 


754 


2 


601 


2 


906 





437 


±0 


042 


2 


752 


2 


601 


2 


943 





329 


±0 


043 


Q2126-158 


3 


059 


2 


906 


3 


211 





310 


±0 


030 


3 


078 


2 


944 


3 


211 





271 


±0 


030 


Q1422+230 


3 


062 


2 


898 


3 


225 





417 


±0 


037 


3 


058 


2 


898 


3 


218 





423 


±0 


034 


Q1422+230 


3 


389 


3 


225 


3 


552 





512 


±0 


045 


3 


382 


3 


218 


3 


552 





496 


±0 


048 


Q0055-269 


3 


092 


2 


928 


3 


257 





381 


±0 


037 


3 


088 


2 


928 


3 


249 





366 


±0 


039 


Q0055-269 


3 


421 


3 


257 


3 


586 





426 


±0 


033 


3 


411 


3 


249 


3 


586 





445 


±0 


041 


Q0000-262 


3 


489 


3 


312 


3 


667 


1 


107 


±0 


086 


3 


708 


3 


607 


3 


808 





705 


±0 


074 


Q0000-262 


3 


845 


3 


667 


4 


023 





758 


±0 


049 


3 


912 


3 


808 


4 


023 





811 


±0 


065 


Q1055+461 


3 


498 


3 


320 


3 


676 





595 


±0 


041 


3 


517 


3 


338 


3 


695 





644 


±0 


056 


Q1055+461 


3 


854 


3 


676 


4 


033 





774 


±0 


050 


3 


862 


3 


695 


4 


033 





843 


±0 


075 


Q2237-061 


3 


880 


3 


690 


4 


070 





828 


±0 


059 


3 


862 


3 


696 


4 


032 





839 


±0 


069 


Q2237-061 


4 


261 


4 


070 


4 


451 





890 


±0 


063 


4 


287 


4 


032 


4 


448 





827 


±0 


060 



B. THE OBSERVED EFFECTIVE OPTICAL DEPTH 



